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Abstract: Today, complex networks have attracted increasing attention from 
various fields of science and engineering. It has been demonstrated that many 
complex networks display various synchronization phenomena. In this paper, 
we introduce a time-varying complex dynamical network model. We then 
further investigate its synchronization phenomenon and prove several network 
synchronization theorems. Especially, we show that synchronization of such 
a time-varying dynamical network is completely determined by the inner- 
coupling matrix, and the eigenvalues and the corresponding eigenvectors of 
the coupling configuration matrix of the network. 
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I. INTRODUCTION 

>- '. A complex network is a large set of interconnected nodes, in which a node is a fundamental unit, 

•i-H . 

^ . which can have different meanings in different situations, such as chemical substrates, micropro- 
cessors or computers, schools or companies, papers, webs, even people, and so on [1-19]. Examples 
of complex networks include the Internet, the World Wide Web, food webs, electric power grids, 
cellular and metabolic networks, etc. [1,2]. These large-scale complex networks often display 
better cooperative or synchronous behaviors among their constituents. 

Traditionally, complex networks were studied by graph theory, where a complex network is 
described by a random graph, for which the basic theory was introduced by Erdos and Renyi [16]. 
Recently, Watts and Strogatz (WS) [4] introduced the concept of small-world networks to describe 
a transition from a regular lattice to a random graph. The WS networks exhibit a high degree of 
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clustering as in the regular networks and a small average distance between two nodes as in the 
random networks. Moreover, the random graph model and the WS model are both exponentially 
growing networks, which are homogeneous in nature. However, according to Barabasi and Albert 
[3], empirical results show that many large-scale complex networks are scale- free, such as the 
Internet, the WWW, and metabolic networks, among others. Notably, a scale- free network is 
inhomogeneous in nature; that is, most nodes have very few connections but a small number of 
particular nodes have many connections. 

In order to better understand the dynamical behaviors of various complex networks, one may 
extend the existing network models from static to dynamic by introducing dynamical elements 
into the network nodes. For the resulting dynamical networks, one significant and interesting 
phenomenon is the synchrony of all dynamical nodes [5-12,15,17-21]. In fact, synchronization 
is a basic motion in nature that has been studied for a long time [5-12,17-21]. Recently, the 
synchronization of chaotic coupled networks has become a focal point in the study of nonlinear 
dynamics [5-12,15,17-21]. We observe that most existing work have been devoted to completely 
regular networks, including lattices, fully-connected graphs, grids, and chains, etc. [5-12,17,21], 
where examples include discrete-time coupled map lattices (CML) [5] and continuous-time cellular 
neural networks (CNN) [13]. A common characteristic of these kinds of regular networks is that 
their complexity originates from the nonlinear dynamics of the nodes, but not from the network 
structures. 

Recently, Wang and Chen [15] introduced a simple uniform scale-free dynamical network model. 
Although their model reflects the complexity from the network structure, it is only a simple uni- 
form dynamical network with the same coupling strength for all connections. In reality, complex 
networks are more likely to have different coupling strengths for different connections. More- 
over, real-world complex networks, particularly biological networks, are time-varying networks. 
Therefore, in this paper, we attempt to introduce a more general time- varying dynamical network 
model, and further investigate its synchronization properties. Based on this model, we prove sev- 
eral network synchronization theorems, showing that synchronization of a time-varying dynamical 
network is completely determined by its inner-coupling matrix (the connecting matrix between 
nodes), and the eigenvalues and the corresponding eigenvectors of its coupling configuration ma- 
trix. On the other hand, many real-world complex networks are directed, for instance the World 
Wide Web. Our new time-varying dynamical network model can well characterize both directed 
and undirected dynamical networks. 

The rest of this paper is organized as follows. A new time- varying dynamical network model is 
introduced in Section II. Section III establishes several network synchronization theorems, along 
with an example for illustration. Conclusions are finally given in Section IV. 
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II. A TIME- VARYING COMPLEX DYNAMICAL 

NETWORK MODEL 



In this section, we introduce a general time- varying complex dynamical network model and several 
mathematical preliminaries. 

A. A General Time- Varying Dynamical Network Model 

Consider a dynamical network consisting of JV linearly and diffusively coupled identical nodes, 
with each node being an n— dimensional dynamical system. The proposed general time-varying 
dynamical network is described by 

JV 

x, = f(xt) + ]T cy(t) A(t) (xj - Xi), % = 1, 2, • • • , JV , (1) 

3=1 

where Xj = (xn, x i2 , ■ • • , x in ) T G R n is the state variable of node i, A(t) = (dki(t)) nxn G R nxn 
is the inner-coupling matrix of the network at time t, C(t) = (<Hj(t)) NxN is the coupling config- 
uration matrix representing the coupling strength and the topological structure of the network at 
time t, in which Cij(t) is defined as follows: if there is a connection from node i to node j (j ^ i) 
at time t, then Cij(t) ^ 0; otherwise, Cij(t) = (j ^ «), and the diagonal elements of matrix C(t) 
are defined by 

N 

cu(t) = i = 1, 2, iV. (2) 

The above time-varying network $IJ) can be rewritten as 

x, = f (xi) + £ A(t) x„ i = 1, 2, ■ ■ ■ , N . (3) 

Suppose that network Q is connected in the sense that there are no isolate clusters, that is, C{t) 
is irreducible. We assume that with their associated inner coupling matrix the oscillators do not 
display Turing [20] bifurcations or other instabilities at larger coupling values and, hence, the 
results in this paper hold only for that class of oscillators and coupling matrices. 

When A(t), C(t) are constant matrices, network (JSJ) becomes a time-invariant dynamical net- 
work: 

N 

^ = f(xj) + J2 c v Ax j> i = 1, 2, JV. (4) 

i=i 

Obviously, the simple uniform dynamical network of Wang and Chen [15], 

Xj = f(xj) + c ^2 Cij Axj, % = 1, 2, • ■ • , JV, (5) 

3=1 

where C is a — 1 matrix and A is a — 1 diagonal matrix, is a special case of network (fl|). 
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B. Mathematical Preliminaries 

Lemma 1: Suppose that an irreducible constant matrix C = (cij) NxN satisfies condition ©. 
Then: 

(i) is an eigenvalue of matrix C, associated with eigenvector (1, 1, 1) T . 

(ii) If Cij > for all 1 < i, j < N (j ^ i), then the real parts of all eigenvalues of C are less 
than or equal to 0, and all possible eigenvalues with zero real part are the real eigenvalue 0. 
In fact, is its eigenvalue of multiplicity 1. 

The proof is omitted here since it can easily be deduced from the Gerschgorin's circle theorem 
and the Perron- Frobenius theorem [21,22]. 

Definition 1: Let Xj(t , X ) (i = 1, 2, • • • , N) be a solution of the nonautonomous dynamical 
network 

±i = f (xj) + gi(t, x x , x 2 , ■ ■ ■ , xjv), i = 1, 2, • • • , JV, (6) 

where X = ((x°) T , ■ • • , (x^) T ) r G R nJV , f : D -> R n and gi : D x ■ ■ • x D -> R n (z = 
1, 2, • • • , N) are continuously differentiable with D C R n , and gi(t, x, x, • • • , x) = for all t. If 
there is a nonempty open subset E C D, with x- G E (i = 1, 2, ■ ■ • , N), such that Xj(£ , X ) G D 
for all t > 0, i = 1, 2, •• JV, and 

lim \\xi(t, Xo) -s(t, xo)|L = for 1 < i < , (7) 

t — >oc 

where s(t , x ) is a solution of the system x = f (x) with x G D, then the dynamical network (jUJ) 
is said to realize synchronization and E X ■ ■ ■ X E is called the region of synchrony for network (JBJ). 
Moreover, X(t,X ) = (xf(£,X ), x^(t,X ), x^(t,X )l is called the synchronous solution 
of network ©, if x;(t, X ) = Xj(t, X ) for all t > and 1 < i, j < N. 

Definition 2: [23] Let T = {s(t) \ < t < T} denote the set of T— periodic solutions of system 
x = f(x) in R™. A T— periodic solution s(t) is said to be orbitally stable, if for each e > 
there exists a 5 > such that every solution x(£) of x = f(x), whose distance from T is less 
than 5 at t — 0, will remain within a distance less than e from T for all t > 0. Such an s(t) 
is said to be orbitally asymptotically stable if, in addition, the distance of x(i) from T tends to 
zero as t — > oo. Moreover, if there exist positive constants a, /? and a real constant h such that 
||x(t — h) — s(t)|| < ae"'* for £ > 0, then s(t) is said to be orbitally asymptotically stable with 
an asymptotic phase. 

For network (J3J), diffusive coupling condition (J2J ensures that (xf (£, X ), x^(t, X ), • • • , x^(t, X )) T 
with Xj(t, X ) = s(t , x ) (1 < z < N) is a synchronous solution of network ©, where 
X = ((x'j 1 ) 2 ", ■ • • , (x^) T ) T , x? = ••■ = x° N = x G -D, and s(£, x ) is a solution of sys- 
tem x = f(x). It is noted that s(t , xo) can generate an equilibrium point, a periodic orbit, an 
aperiodic orbit, or a chaotic orbit in the phase space. In the following, we only study the case 
that s(t , x ) cannot generate a chaotic orbit, since we have already investigated the case of chaos 
synchronization in general time- varying dynamical networks in [18]. 

Note that non-chaotic synchronization and chaos synchronization are different in essence. It 
is very important to point out that since we do not assume the stability of s(t , xq) for system 
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x = f(x), the stability of the synchronous solution S(£) = (s T (t, x ), s T (t, x ), ■ ■■ ,s T (t, x )) T 
of network (jSJ) is equivalent to the stability of the error vector ^(t), ■ • • , Vn(^ about 

its zero solution, where rji(t) = Xj(t) — s(£ , x ) {i = 1, 2, • • • , N). However, it is quite dif- 
ferent from the chaos synchronization case, since a chaotic attractor is an attracting invariant 
set and a chaotic orbit is not stable in Lyapunov sense, so the stability of the chaotic syn- 
chronous state xi(t) = x 2 (t) = ••• = xjv(f) is equivalent to the stability of the transverse 
errors (^(t), vj(t)i *"> Vn(^ °f the synchronous manifold of network (JHJ), where r]i(t) = 
Xj(t) — s(t , x ) {i = 2, 3, • ■ ■ , N) and x^t) = s(t , x ) is the reference direction of the syn- 
chronous manifold [18,19]. 



III. SEVERAL NETWORK SYNCHRONIZATION 

CRITERIA 

In this section, we establish several network synchronization criteria for both the time-varying 
network (JHJ) and the time- invariant network 



A. Periodic Orbits Synchronization of Time-Invariant Dynamical Net- 
works 

Here, we consider periodic orbits synchronization of the time-invariant network (0J). Assume that 
s(t , x ) is a nonconstant T— periodic solution of an individual node x = f(x), denoted s(t). 
Now, let 

rjiit) = Xi(t) - s(t), i = 1, 2, N, (8) 
and substitute (JBJ) into network (jlj) to get 



N 



Ut) = f(s(t) + r)i(t)) - f(s(t)) + £cy Atfe(t) , 1 < % < N. 



Also, denote 

/ *a(t) \ / Vi(t) \ 

x(t) = ; , v(t) = \ 

V *N(t) J V VN(t) j 

Then, Q can be rewritten as 

f)(t) = F(t, fj(t)) 
and the Jacobian of F(t, fj) at fj — is 

/ Df(s(t)) + c n A 



S(t) 



/ s(t) \ 
V s(t) y 



G R 



nAT 



L>F(t, 0) 



Pai A 

Cjvi A 



Cl2 A 

Df(s(t)) + c 22 A 

CAT3 A 



civ A 
c 2 jv A 

Df(s(t)) +c NN A J 



(9) 



(10) 
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Obviously, S(t) = (s T (t), s T (t), ■ ■ • , s T (t)j is a synchronous periodic solution of network (j3J). 
It is interesting and also important to find out whether or not this synchronous periodic solution 
S(t) is also orbitally asymptotically stable in network (j3J) if s(i) is orbitally asymptotically stable 
in an individual node described by x = f(x). The following theorem gives a positive answer to 
this question. 

Definition 3: Suppose that s(t) is a periodic solution of system x = f (x). Let 71 = 1, 72, • • • , 7 n 
be the Floquet multipliers of the variational equation of s(£), y = A(t) y, where A(t) = Dt(s(t)) 
is the Jacobian of f evaluated at s(t). Then the periodic solution s(t) is said to be a hyperbolic 
periodic solution if 7^ 1 for 2 < j < n. Moreover, S(t) is said to be a hyperbolic synchronous 
periodic solution of network (j3J) if all Floquet multipliers of the variational equation of S(t) have 
absolute values less than 1 except one multiplier which equals 1. 

Theorem 1: Suppose that s(t) is a hyperbolic periodic solution of an individual node x = f (x), 
and is orbitally asymptotically stable with an asymptotic phase. Suppose also that the coupling 
configuration matrix C = (cij) NxN can be diagonalized. Then, S(t) is a hyperbolic synchronous 
periodic solution of network Q, and is orbitally asymptotically stable with an asymptotic phase, 
if and only if the linear time-varying systems 

w = [Df (s(*)) + \ k A] w , k = 2, ■ ■ ■ , N , (11) 

are asymptotic stable about their zero solutions. 

Proof. It is well known that the stability of a periodic solution of an autonomous system is 
completely determined by its Floquet multipliers [23]. We linearize network (JH) at X(£) = S(t) 
and get 

Y(t) = DF(t, 0) Y(t) , (12) 

where DF(t, 0) is defined in dTUJ) and Y(t) = (yj(t), yJ?(t)) T G K nN . According to the 
stability theory of periodic orbits [23], S(t) is a hyperbolic synchronous periodic solution of network 
(@J), and is orbitally asymptotically stable with an asymptotic phase, if and only if all Floquet 
multipliers of system (fT2*|) have absolute values less than 1 except one multiplier which equals 1. 
From (|12jl . we have 

N 

y t (t) = Df(s(t))yi(t) + E fkj Ay j(t) 

j'=i 

= Df(s(t)) y<(t) + A ( yi (t), • • • , y N (t)) (c a , • • • , c iN ) T , 1 = 1, 2, • • • , N , 

where Y(t) = (y x (t), • • • , y N (t)) G K nxN . That is, 

Y(t) = Df(s(t))Y(t) + AY(t)C T . (13) 

Obviously, systems f!12l) and IjlHj) have the same Floquet multipliers, since they are different only 
in form. 

Since C can be diagonalized, namely, it has all real eigenvalues of multiplicity one, there exists 
an nonsingular matrix, <I> = (0i, • • • , 4>n), such that 

C T $ = $A, 
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where A = diag{Ai, A2, • ■ • , Ajv}, and Ai, • • • , Aat are the eigenvalues of C. Let 

Y (t) = V(t) $ 1 . 

From (JT3J), the matrix vector V(t) = (vi(t), • • • , Vjv(t)) € R nxAr satisfies the following equation: 

V(t) = Df(s(*))V(t) + AV(t)A, 

namely, 

v fc (t) = [Df (s(*)) + A, A] v fc (t) , fc = 1, 2, • • • , A . (14) 

Thus, we have changed the stability problem of the synchronous periodic solution S(t) into the 
stability problem of iV independent n— dimensional linear systems ()14j) . which has the same form 
as (fTTTl . 

It follows from Lemma 1 that one eigenvalue of C satisfies Ai = and the corresponding 
linear system in (fT4*j) is 

vx(£) = Df(s(t))vi(t). (15) 

Obviously, system (fl3|) is the corresponding linearized system of an individual node x = f(x) at 
x = s(t). When k > 1, systems (fT4*|) become systems (JTTJ). 

Let 7ij (i = 1, 2, • ■ • , A; j = 1, 2, - • • , n) be the Floquet multipliers of the A independent 
n— dimensional linear time- varying systems (j!4)l . Obviously, the Floquet multipliers of system 
((H) (or dU)) are also given by 7^ (z = 1, 2, ■ ■ ■ , N; j = 1, 2, • • • , n). 

According to the stability theory of periodic orbits [23], s(t) is a hyperbolic periodic solution of 
an individual node x = f(x), and is orbitally asymptotically stable with an asymptotic phase, if 
and only if all Floquet multipliers of (|T3j) satisfy \ < 1 for j = 2, • • • , n, except one multiplier 
711 which equals 1. Moreover, the linear systems (fTTjl are asymptotically stable if and only if 
all Floquet multipliers of systems (|TT|l satisfy |7y| < 1 for % = 2, • • • , N and j = 1, 2, • • • , n. 
Therefore, S(t) is a hyperbolic synchronous periodic solution of dynamical network Q, and is 
orbitally asymptotically stable with an asymptotic phase, if and only if all Floquet multipliers of 
the linear systems (fT4"|) satisfy | 'y ij I < 1 for j = 2, • • • , n and < 1 for i = 2, • • • , N; j = 
1, 2, ■ ■ • , n, except 711 = 1. 

The proof is thus completed. 

Note that the result of Theorem 1 can not be extended to the time- varying dynamical network 
Q, which is a nonautonomous system. Network (j3J) is discussed in the next subsection. 

B. Stable Orbits Synchronization of Time- Varying Dynamical Networks 

Here, assume that s(t , x ) is an exponentially stable solution of x = f (x), simply denoted as s(£). 
Let fi (A) be the maximum eigenvalue of \ (A T + A) . We provide two network synchronization 
theorems for the time-varying dynamical network (j^J) in this subsection. 
Substituting (JHJ) into the time-varying network (jSJ) yields 

JV 

fu(t) = f (s(t) + - f (s(t)) + ]T ^{t) A(t) rijit) , 1 < % < N . (16) 
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Similar to (jlOj). system ()16|) can be reformulated as 

fj(t) = F(t, fj(t)) 
and the Jacobian of F(t, fj) at fj = is 

/ Df(s(t)) + c n (t)A(t) c 12 (*)A(t) 



DF(t,0) 



Q2l(*)A(t) 
CATl(t)A(t) 



Df(s(t))+C2 2 (t)A(t) 

c N2 (t)A(t) 



(17) 



c 1Ar (t)A(£) \ 
c 27V (t)A(t) 

Df (s(t)) + CAW (i)A(i) y 



Theorem 2: Suppose that F : 0^ R nAf is continuously differentiable on = {fj G R nAf |||f/|| 2 < 
r}. The synchronous solution S(t) is uniformly exponentially stable in the dynamical network (JHJ) 
if there exist two symmetric positive definite matrices, P, Q G pj lNxnN ; such that 



P [DF(t, 0)J + (DF(t, 0)J P < -Q < -cil, 

where c\ > 0, and 

(r(t, y ) - r(t, s(t))) T p + p (r(t, y) - r(t, s(*))) < c 2 i< Cl i, 

in which T(t, y(t)) = diag{Df ( yi (t)), • -• , J»(y*(t))}, y(t) = (yf (*), yJC*). ■ • • , y T N 
y-S(t) G fi. 

Proof. Since F is continuously differentiable on Cl, from ([16 )1 - ([17 )1 . we have 

/ f(m(t) + s(*))-f(s(t)) \ / c u (t)A(t) •■■ c 1JV (t)A(t) \ 



and 



f( VN (t) + b(*)) - f(s(t)) y \ CJV i(t) a( 

Df( yi (t,^(t)))^(t) \ / c u (t)A(t) 



7j(t) 



• • c^(t)A(*) y 

c liV (t)A(t) \ 



c NN (t)A(t) ) 



Df(y N (t, VN(t))) VN(t) J \c m (t)A(t) ■ 
Df( yi (t, m (t))) - Df(s(t)) 

'•. 

Df(y N (t, rj N (t))) - Df(s(t)) y 



»7(t) 



+ £>F(t, 0) 77(t) 
= r(t, y(t)) - T(t, S(t)) + DF(t,0)} fj(t) . 

Define a Lyapunov function, using the vector fj(t), by 

= m T p m , 

and differentiate V(t) with respect to time £ to get 
V(t) = rj(t) T Pfj(t) + fj(t) T Pfj(t) 



fj(t) 



T 



fj{t) 



DF(t, 0)j P + P (DF(t, 0) 

+ n{t) T [(r(t, y ) - r(t, s(t))) T p + p (r(t, y ) - r(t, s(t)))J r,(t) 

< -fj(t) T Qfj(t) + c 2 fj(t) T Ifj(t) 

< (Pa - ci) fj{t) T fj{t) < 0. 



8 



From the Lyapunov stability theory, the synchronization errors fj(t) will uniformly exponentially 
converge to zero. That is, the synchronous state S(t) is uniformly exponentially stable for the 
dynamical network Q. The proof is thus completed. 

Next, we consider a class of time- varying dynamical networks whose configuration matrix 
satisfies a special property. 

Assumption 1: Let Ai(t), \2(t), Aat(£) be the eigenvalues of C(£). 3t > 0, for any 
Xi(t) (1 < % <N), either X l (t) ^ for all t > t , or X^t) = for all t > t . 

For any fixed ti > t > C(ti) is a constant matrix. According to Lemma 1, there exists a unique 
Xi(t) satisfying Xi(t\) = and Aj(ti) < for j ^ i. From Assumption 1, the eigenvalues of C(i) 
satisfy Xi{t) = for all t > to and Xj(t) < for all t > to, j ^ i. Let Xi(t) be Xi(t). 

Theorem 3: Let x = s(t) be an exponentially stable solution of nonlinear system x = f(x), 
where f : Q — > R" is continuously differentiable, Q = {x G R n |||x — s(£)|| 2 < r}. Suppose 
that the Jacobian DF(t, ff) is bounded and Lipschitz on = {fj G R n7V \ \\f)\\2 < r}, uniformly 
in t. Suppose also that Assumption 1 holds and there exists a real matrix, 3>(t), nonsingu- 
lar for all t, such that (C(t)) T = diag{Ai(t), A 2 (t), X N (t)} and = 

diag{/3i(f), ffait), ■ ■ ■ , /3;v(t)}. Then, the synchronous solution S(t) is exponentially stable in dy- 
namical network (jSJ) if and only if the linear time-varying systems 

w = [Df (s(t)) + A fc (t) A(t) - /3 fc (t) I„] w , k = 2, • ■ ■ , N , (18) 

are exponentially stable about their zero solutions. 
Proof. Let 

x(t) = z(t) + s(t) , (19) 
and substitute ()19)) into x = f(x), so as to obtain 

z = Df(s(t))z + h(t, z) , (20) 

where h(t, z) = f(z + s(t)) - f(s(t)) - Df(s(t))z. 

Since DF(t, ff) is bounded and Lipschitz on 0, uniformly in t, -Df(x) is also bounded and 
Lipschitz on f2, uniformly in t. That is, the Jacobian of system ()20|). i.e., 8f( - z ^ z s ^' > = d g^+s(t)) = 
-Df(x), is bounded and Lipschitz on Q, uniformly in t. According to the Lyapunov converse 
theorem [24], the origin is an exponentially stable equilibrium point for the nonlinear system 
f!20[) if and only if it is an exponentially stable equilibrium point for the corresponding linear 
time-varying system 

y = Df(s(t))y. (21) 

Since network (j^J can be rewritten in the form of system ()17)) . we linearize network (j3J) at 
X(t) = S(t) and get 

Y = DF (t, 0) Y, (22) 

where Y(t) = (yf(t), • • ■ , y T N (t)) T G R nN . 

Since DF(t, ff) is bounded and Lipschitz on Q, uniformly in t, according to the Lyapunov 
converse theorem [24], the origin is an exponentially stable equilibrium point for the nonlinear 
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system ()17|) if and only if it is an exponentially stable equilibrium point for the corresponding 
linear time- varying system (|22|) . Note that the exponential stability of the zero solution of system 
(122(1 means that Y(£) — > exponentially. Obviously, Y(t) — ■> exponentially is equivalent to 
Y(t) -> exponentially, where Y(t) = ( yi (t), y 2 (t), • ■ • , yjv (t)) G R nxiV . 

Let Y(t) = V(£) <fr _1 (t), where <fr(t) is defined in the statements of Theorem 3. Similar to 
the proof of Theorem 1, we have 

v fc = [Df (s(£)) + X k (t) A(t) - /3 k (t) I n ] v fe , k = 1, 2, • • • , N . (23) 

Thus, we have changed the exponential stability problem of the synchronous solution S(t) into 
the exponential stability problem for iV independent n— dimensional linear systems (|23|). which 
has the same form as (jl8)l . 

From Assumption 1 and Lemma 1, \\{t) = for all t > to, and so the corresponding system 
(J23|) becomes system (|21|) . From the assumed conditions of Theorem 3, vi(t) — > exponentially 
Moreover, Y(t) — > exponentially is equivalent to Vfc(t) — > exponentially for A; = 2, • ■ ■ , N. 
Therefore, the synchronous solution S(t) is exponentially stable for network © if and only if the 
linear time- varying systems (jl8|) are exponentially stable about their zero solutions. The proof is 
thus completed. 

Remarks: Theorem 3 shows that synchronization of the time- varying dynamical network © is 
completely determined by its inner-coupling matrix A(i), and the eigenvalues A&(t) (2 < k < N) 
and the corresponding eigenvectors ^ (t) (fi^it) are functions of , 2 < k < N) of the 
coupling configuration matrix C(t). However, the synchronization of the time-invariant dynamical 
network (0J) is completely determined only by its inner-coupling matrix A and the eigenvalues of 
the coupling configuration matrix C [14,15]. 



C. A Simulation Example 



In this subsection, we illustrate Theorem 3 by using a 3-dimensional exponentially stable system 
as a node in network (JHJ). For simplicity, we only consider a three- node network. Each individual 
node is described by X\ — — Xi, X2 — — 2 X2-, A3 — — 3X3, which is exponentially stable at 
s(£) = 0, and its Jacobian is -Df(x) = diag{ — 1, —2, —3}. 

Assume that the inner-coupling matrix is A(t) = diagjl + e~ 2t , 1 + e~ 3t , 1 + e~ 4 }. and the 
coupling configuration matrix is 



C(t) 



2e 2 - e - 1 



/ cn(t) c 12 (t) c l3 (t) \ 
C21O) c 22 (t) c 23 (t) 

V C3i(t) C 32 {t) C 33 (t) J 



where cn(t) 



l)th(t) + e arctan(t), ci 2 (t) 



e) th(t) — 2e arctan(t), c^t) 



(e - e 2 )th(t) + e arctan(t), c 21 {t) = 2(e 2 - l)th(t) + e 2 arctan(t), c 22 (t) = 2(1 - e) th(t) - 
2e 2 arctan(t), c 23 (t) = 2(e — e 2 ) th(t) + e 2 arctan(t), c 3 i(t) = 3(e 2 — 1) th(t) + arctan(t), c 32 (t) = 
3(1 - e)th(t) - 2arctan(t), c 33 (t) = 3(e - e 2 ) th(t) + arctan(t), with th(t) = 4x^r- 
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It is easy to verify that there exists a nonsingular real matrix, 

2 2 (1 - e 2 )e* - e 1+t6 »® \ 
1 - 3e (e - lie* 2e 1+sin W 

,2 



*(*) 



1 



2e 2 



e - 1 



\2e - e 2 



(e 2 - e)e* -e 1 + sin « 



such that (C(t)) T $(t) = diag{0, -th(t), -arctan(t)} and *(*) = diag{0, -1, 

— cos(t)}. Obviously, the conditions of Theorem 3 are satisfied. Therefore, the zero solution of 
network (J3J) in this example is exponentially stable if and only if the linear time-varying systems 



w = [£>f(s(t)) + X k (t)A(t) - (3 k (t)I 3 )w, k = 2, 3, 



(24) 



are exponentially stable about their zero solutions. 

When k = 2, wehaveA 2 (t) = -th(t), f3 2 (t) = -1. Then, we have, for all t > 1, n [Df(s(t)) + 
A fc (i)A(i) -^ fc (t)Is] = - (1 + e~ 2 *) th(t) < -th(l) < 0, showing that the decoupled system pEjl 
with A; = 2 is exponentially stable about its zero solution. 

When & = 3, we have \%(t) = — arctan(t), (3^{t) = — cos(t). Similarly, we can verify that 
system (|2*4*|) with A; = 3 is also exponentially stable about its zero solution. 

Therefore, from Theorem 3, the synchronous solution S(t) = of network (J3J) is exponentially 
stable. 

Note that, in this example, C(t) is not a constant matrix, as assumed in other investigations 
[7-12,15,17,21]. 



IV. CONCLUSIONS 

Over the last few years, the discovery of small-world and scale-free properties of many complex 
dynamical networks has seen great advances, while these networks all have constant and uniform 
connections. In this paper, we have further introduced a time-varying dynamical network model, 
and investigated its synchronization criteria. Several network synchronization theorems have been 
established for this model. Especially, we have shown that: synchronization of a time-varying 
dynamical network is completely determined by its inner-coupling matrix, and the eigenvalues 
and the corresponding eigenvectors of the coupling configuration matrix, differing from that for 
time-invariant complex networks studied elsewhere before. 

The proposed time-varying complex network model provides a mathematical description for 
further investigating the dynamical behaviors and topological structures of many real-world com- 
plex dynamical networks. We foresee that this model will be very useful for the current intensive 
studies of general complex dynamical networks. 
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